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1. Introduction 

The standard approach used in the evaluation of form factors in lattice QCD is to compute a 
three-point function. More detailed information on hadron structure can be extracted from four- 
point correlators. The quark distribution inside the hadron and hadron deformation are just two 
such important aspects that can be studied using these correlators. The equal time density-density 
correlator provides a gauge invariant definition of the hadron "wave function" but originally could 
only be evaluated approximately This is because four-point functions are harder to compute 
than two- and three-point functions, requiring the all to all propagator. The usual way to estimate 
the all to all propagator is by employing stochastic techniques [^. In Ref we used Z(2) noise 
combined with dilution to compute the all to all propagators and obtained results both for mesons 
and baryons [||, 0]. In this work we generalize the so called one-end trick [||], originally devised as a 
method to calculate two-point functions, to evaluate four-point functions. We demonstrate that this 
approach yields more accurate results by evaluating the density-density correlator for the pion and 
p -meson and comparing the results to those obtained using standard stochastic techniques [^, ^. 
Furthermore, we extract the pion form factor obtaining results that have comparable errors as those 
obtained when one uses the one-end trick to compute the pion form factor using the three-point 
function [^. An advantage of using four-point functions is that we only need one set of stochastic 
propagators to extract the form factor for any momentum transfer unlike using three-point functions 
where a new set is needed for every momentum. We also show how to generalize our method to 
other mesons and give preliminary results on Gi, one of the three form factors of the p— meson. 

2. Four-point functions 

Hadron four-point functions are given by 

{X2,ti,t2) = J d\id^x{h{x,t)\fJ{x2+xi,t2)ja'{xi,ti)\h{xo,to)) (2.1) 

where is the normal ordered electromagnetic operator : qfjaqf '■ with / being a flavor index, 
while \h) denotes any hadronic state. The two integrations ensure zero momentum of the hadronic 
state; integrating over x\ sets the momentum of the source equal to that of the sink and integrating 
over X sets both to zero. Thus to compute the four-point function on the lattice, the all to all 
propagator from all sites x\ to x is needed. 

It is well known that an estimate for the all to all propagator can be obtained using stochastic 
techniques In brief, one inverts for a set of Nr noise vectors obeying (^^fi{x)£,v\yyj = 5{x — 

y)Sa,b^li.v and {^^{x)'j^ = and estimates the all to all propagator by averaging the product of the 
solution vectors with the noise vectors over the stochastic ensemble. Namely the quark propagator 
Gv4i{x,y) — > (^py{x)^^^ (y)'^ , where ^ is a noise vector and the solution vector. One, therefore, 
replaces every occurrence of G with the product between ^ and <p thereby obtaining the stochastic 
estimate for the four-point function. More details and results for hadron wave functions and the 
pion form factor obtained using this method can be found in Refs. ||3|, ^], where it is shown that 
with sufficient number of noise vectors and dilution one can obtain a reasonable signal [Q]. We 
shall refer to this approach as Method I. Here we show how one can reduce stochastic noise by 
implementing the one-end trick for the computation of meson four-point functions. We shall 
refer to this new approach as Method II. 
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3. Description of Method II 



The one-end trick was originally devised for the precise calculation of pion two-point func- 
tions. In its original form, one combines appropriately solution vectors so that an automatic sum- 
mation over the source coordinate arises. Thus the number of stochastic inversions needed is re- 
duced to a few inversions, thereby suppressing stochastic noise. More explicitly, expanding the dot 
product between two solution vectors yields the pion two-point correlator summed over the source 
coordinate: 



x,yo,xo 



where we assume that the noise vectors are localized on a certain time slice tQ. Taking the average 
of the noise vectors over the stochastic ensemble yields 5— functions by definition. Thus we obtain 



^ Gfy{x,t;xo,to) Gfy,{x,t;yo,to)5aa'5vv'5{xo-yo) = J^Tr \G{x,t;xo,to)\ 



(3.2) 



x,yo,xo 



x^o 



In the case of the pion Eq. (3.2) arises automatically since one combines the backward going prop- 
agator with the 75 pairs that appear in the pion interpolating operator. For a general interpolating 
field, say qiTqj, where Y is any product of y-matrices, the noise vectors must be diluted in spin 
i.e. t,^{xQ,tQ)r = ^''{xo,to)5r^, r = 1,...,4. This imposes that the number of noise vectors is in 
multiplets of four. In this case the solution vectors give (/)^(x,f;fo)r = Ljfo G"^^{x,t;xo,to)^''{xo,to), 
where the to argument appearing in (p is to remind us that the noise vector is localized on the time 
slice to- The combination given by 



£(/>;(x,f;fo)v(r75)vr0r(^,f;fo)r(75n-cM, 



(3.3) 



yields the two-point function of the meson summed over both sink and source coordinates. The 
downside of this method is that, due to the automatic summation over both the source and sink 
spatial coordinates, one cannot compute two-point functions for arbitrary momenta using a given 
set of noise vectors. To utilize the one-end trick and extract the two-point correlator at finite mo- 
mentum one must multiply the noise vectors with an appropriate phase prior to the inversion. Thus 
one needs Nr inversions for every momentum vector. 

The application of this method to four-point functions is appealing since by definition (see 
Eq. ( |2.1[ )) one is interested in the case where the initial and final states are at rest. We show here 
how the one-end trick can be implemented for the case of mesons. The four-point function that we 
consider is shown schematically in Fig. |l|. At the propagator level we have 



G{?{x2;to,h,t2,t)= £ {Xr{x,t)\f{x2+xi,t2)f{xuh)\xrixQ,tp)) (3.4) 

Xl ,X,XQ 

Tr[ys yaG{xi ,h;xo, to)Yy5 G^ (f 2 + -^i , ^2 ; -^o , ^o) 75 yaG{x2 +xi,t2;x, t)Tys G^ {x^,ti;x,t)] 



Xl ,x,xo 



where Y = yoF^yo and Xt{x) is the interpolating field of the meson that takes the general form 
qi{x)Yqj{x) for / / j. Here we explicitly add a sum over the source coordinate xq and we fix the 
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time slice of the source, to, and the sink, t. The time 
slices, ti and t2, where the currents are inserted, on the 
other hand, are free to take any value between to and 
t. Thus one needs two sets of stochastic inversions, 
one set with the noise vectors localized on the time 
slice, to and one on the time slice, t. One then finds 
an appropriate combination of solution vectors such 
that the summation over source and sink coordinates 
is carried out automatically. The combination: 

Y^Tr [Y5'Yc!S{T;xuti;x2+xi,t2;to)Y5YaSir;x2+xi,t2;xi,ti;t)] (3.5) 
where S'^y{r;x2+xi,t2;xi,ti;t) = Lr0^'(^2 +-?i,^2;0f (r75)ric0v^(-^i,^i;Oic achieves this. 

Throughout this work we use two degenerate flavors 
of dynamical Wilson quarks [§,0. In all computa- 
tions we employ Gaussian smearing combined with 
hypercubic (HYP) smearing of the gauge links that 
enter the Gaussian smearing function. The parame- 
ters of the Gaussian smearing are adjusted to ensure 
minimal time evolution for filtering the meson ground 
state. The parameters of our calculation are summa- 
rized in Table |l|. 

The p -meson charge distribution is obtained using the equal time density-density correlator 
given by 

Gp{xi,ti) = J d^X2d^x{p{x,t)\f^{x2+xuti)j^{x2,ti)\p{xo,to)) . (4.1) 

We test the new technique by comparing results for Gp (^i , ti ) using Methods I and II. In the large 
?i and t — ti limit when the p state dominates, G^°{xi,ti), normalized over the spatial volume, 
becomes time independent and it is denoted by Cp{xi). In the non - relativistic limit, this four- 
point function reduces to the wave function squared. The ingredients needed in Method I are the 
point to all propagator from the source and two all to all propagators at time slices ti and t, both of 
which are kept fixed. For this computation we use six sets of noise vectors diluted for each spin, 
color and even-odd sites i.e. we need 24 x 6 = 144 inversions to obtain each stochastic propagator. 
This means a total of 144 x 2 + 12 = 300 inversions are required for each gauge configuration 
For Method II, on the other hand, we used eight sets of spin diluted noise vectors at the source 
and sink thus a total of 64 inversions for each configuration. In Fig. ^ we show a comparison 
between the results obtained using Method I and II. What is plotted is the projection of the density- 
density correlator, Cp(xi), along the spin axis taken to be the z-axis and perpendicular to it. The 
interpolating field used for the p -meson is UY3d. As can be seen the statistical errors obtained when 



j3 =5.6, fl-' =2.56(10) GeV 



# Confs. 


K 


am. 


r 


m;r/mp 


24^ X 40 [8] 




185 
150 


0.1575 
0.1580 


0.270(3) 
0.199(3) 


0.69 
0.56 


24^ X 32 [9] 




200 


0.15825 


0.150 


D) 


0.45 



Table 1: The simulation parameters used in 
our computations. 



f (X2+Xi,?2) 




Figure 1: The four-point function for 
mesons. 

4. Meson wave functions 
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Figure 2: Projections of the p -meson 
density-density correlator along the spin axis 
and perpendicular to it. Upper for Method 1 
and lower for Method 11 at K = 0. 1580. 



using Method II are almost four times smaller despite 
the fact that we use 144/32 = 4.5 less number of noise 
vectors to estimate the all to all propagator. Therefore 
the improvement gained using the one-end trick is re- 
ally significant, reducing computational time by two 
orders of magnitude. The results obtained in Method 
II clearly reveal an asymmetry in the charge distribu- 
tion of the p -meson, which in Method I was hard to 
see. 

Having demonstrated the effectiveness of Method II 
we use it, in what follows, to study deformation in the 
p -meson as a function of the quark mass and to ex- 
tract the pion and p-meson form factors. In Fig. || we 
show contour plots of the density-density correlator of 
the p-meson, Cp(r), projected onto the y-z plane. As 
can be seen, for all three pion masses, we obtain an 
ellipse that is elongated along the spin axis, showing a 
clear deformation from spherical symmetry. The cor- 
responding contour plots for the pion show no devia- 
tion from the circle as expected. 




Figure 3: Contour plots of the charge distribution of the p -meson projected onto the y-z plane for 
all three K values studied. The dashed circles are to guide the eye. 



5. Pion and p-meson form factors 

Form factors can be accurately extracted using four-point functions by taking the Fourier trans- 
form of G-j^" {x,ti,t2) and allowing large time separations between the current insertions, t2 — ti. 
Therefore the extraction of form factors requires larger temporal extension than the equal time 
density-density correlators. Methods to suppress excited state contributions are therefore of crucial 
importance here. Gaussian smearing combined with HYP smearing achieves ground state domi- 
nance as early as three time slices. Taking the Fourier transform of the pion four-point correlator 
we obtain 
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1^75 1^(0)) 



.2\{E{p)+mMQ^)\ 



.g-E(p){t2-ti)^-m„(t-{t2-ti)-to) 



(5.1) 



%mlE{p) 

where Fji is the pion form factor and is the Euclidean momentum transfer squared. The time 
dependencies and overlaps cancel by dividing with an appropriate combination of two-point func- 
tions: 



Ry,{p;tut2) 



g4 {p;tQ,ti,t2,t)Gy^{p,ti-tQ) 



E{p)+m^ ^ G^^{p^t2-tQ)Gy,^,t-{t2-h)-to) 



(5.2) 



where Gy^{p,t) is the pion two-point function at momentum p. We search for a plateau of 7?^^ {pUiJi) 
by varying the time difference t2 — ti, as shown in Fig. |^ We perform the calculation for two source 
- sink separations, namely {t — to) / a =14 and 16 to check that we have ground state dominance. 
As can be seen, we obtain consistent plateau values for both source - sink separations. 
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Figure 4: Ry^ {p;ti,t2) versus {t2~~t\)/a for K = 
0. 1575. The range used for the fit is shown by the 
length of the lines. 




3.0 



1.0 1.5 2.0 

Q2 (GeV2) 

Figure 5: The pion form factor for three K- 
values. We compare with results using the hy- 
brid approach from [|T^ and twisted mass results 
from [|. 

In Fig. ^ we show the pion form factor for three K"- values compared with recent results obtained 
using three-point functions. Results in the hybrid approach, that uses dynamical staggered sea 
quarks and domain wall valence quarks, are obtained using sequential inversions to compute the 
three-point function [|I0|]. Results with dynamical twisted mass fermions, on the other hand, use 
the one-end trick to compute the three-point function [^. Our results compare very well to those 
obtained in the latter case, which is closest to our approach. Assuming vector meson dominance 
and taking lUp = 0.77 GeV we obtain the curve shown, for reference, in Fig. ^. 

The p -meson has a Coulomb, Gc, a magnetic, Gm and a quadrupole, Gq form factor. They 
can be parametrized in terms of Gi, G2 and G3 as 



Go = Gi-G2+ 1 + 



21 
4M2 



G3, 



Gm — G2, Gc 



2 



One can find combinations between initial and final p polarizations and insertion directions (a) 
that isolate individual form factors and for which decay to a pion is forbidden. Here we shall 
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demonstrate the method by showing preUminary results for Gi. As for the case of the pion form 
factor, Gi, can be extracted using only jo: 

X\0) 



Gi {pj_Ui,t2) 



SM^E{pA 



{E{p^)+MYG\{Qy 



-M{t-(t2-tl)-to)^-E{pj_){t2-tl) 



(5.3) 



1.0!! 
0.8 ■ 
0.6 ■ 
0.4 ■ 
0.2 ■ 
0.0 ■ 



Gi(Q2) 



OTp = 1.002(41) GeV h 
nip = 0.910(38) GeV h 
nip = 0.853(37) GeV i- 



■ 1 

Q2 (GeV^) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
Figure 6: Gi as a function of the momentum transfer 
for three K values. 

6. Conclusions 



where {Q.\Xr,\p{p,s)) = X{p)£k{p,s), 
Ls£kip,s)£l,{p,s) = gkk' - and pi_ 

is a momentum perpendicular to the k direc- 
tion. As in the case of the pion form factor, we 
construct an appropriate ratio and search for 
a plateau in t2 — t\. Results for G\ are shown 
in Fig. m They carry small statistical errors 
demonstrating the applicability of the method 
for the extraction of the p -meson form factors. 
An analysis to extract all three form factors 
and subsequently derive physical quantities is 
in progress. 



We have shown that the one-end trick can be applied to evaluate accurately four-point func- 
tions. Using this approach, the density-density correlators are computed to sufficient accuracy to 
show that the p -meson is deformed. We also obtain accurate results for the pion form factor that 
compare favorably to the accuracy obtained using the one-end trick to compute the three-point 
function. The advantage of using four-point functions is that only one set of inversions is needed 
for all momentum transfers, unlike in the case of three-point functions where one needs new inver- 
sions for each value of the momentum transfer. Preliminary results on the p -meson form factor, Gi , 
carry small statistical errors demonstrating the applicability of the method also in the calculation 
of the form factors of the p -meson. 

References 

[1] C. Alexandrou, Ph. de Forcrand and A. Tsapalis, Phys. Rev. D 66 (2002) 094503. 
[2] C. Michael and J. Peisa, Phys. Rev. D 58 (1998) 034506. 
[3] C. Alexandrou, P Dimopoulos, G. Koutsou and H. Neff, PoS (LAT2005) 030. 
[4] C. Alexandrou, G. Koutsou and H. Neff, PoS (LAT2006) 113. 
[5] C. McNeile and C. Michael, Phys. Rev D 73 (2006) 074506. 
[6] S. Simula et al, PoS (LAT2007) 374. 
[7] J. Foley et al, Comput. Phys. Commun. 172 (2005) 145. 
[8] B. Orth, T. Lippert and K. Schilling, Phys. Rev D 72 (2005) 014503. 
[9] C. Urbach, K. Jansen, A. Shindler and U. Wenger, Comput. Phys. Commun. 174 (2006) 87. 
[10] F Bonnet et al, Phys. Rev D 72 (2005) 054506. 



7 



